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The density-matrix renormalization-group method is ap- 
plied to the one-dimensional Kondo-lattice model with the 
Coulomb interaction between the conduction electrons. The 
spin and charge gaps are calculated as a function of the ex- 
change constant J and the Coulomb interaction U c - It is 
shown that both the spin and charge gaps increase with in- 
creasing J and U c ■ The spin gap vanishes in the limit of J — > 
for any U c with an exponential form, A s oc exp [— l/a(U c ) Jp]. 
The exponent, a(U c ), is determined as a function of U c . The 
charge gap is generally much larger than the spin gap. In 
the limit of J — > 0, the charge gap vanishes as A c = \ J for 
U c — but for a finite U c it tends to a finite value, which is 
the charge gap of the Hubbard model. 

PACS numbers: 75.30.Mb, 71.28.+d, 71.30.+h 



Recently the insulating phase of the Kondo-lattice 
model (KLM) has been extensively studied in relation 
to the Kondo insulators Q. In contrast to the strong- 
coupling limit where various properties are easily under- 
stood from the local bases, the weak-coupling limit has 
not yet been fully understood. Since the KLM with weak 
exchange coupling is an effective model of the symmet- 
ric periodic Anderson model with strong Coulomb inter- 
action, the weak-coupling region of the KLM is impor- 
tant to understanding the strongly correlated insulating 
phase. For the one-dimensional case it has been estab- 
lished that both spin and charge gaps exist for any ex- 
change coupling. The excitation gaps vanish in the limit 
of vanishing exchange coupling [[2|j3|] . Since the gaps are 
tiny in the weak-coupling regime, it is no longer justified 
to neglect Coulomb interaction between the conduction 
electrons. 

In this article we determine the spin and charge gaps 
in the entire region of exchange constant taking account 
of the Coulomb interaction between the conduction elec- 
trons U c . The excitation gaps are obtained precisely 
by employing the density matrix renormalization group 
(DMRG) method. Among other things, we find that U c 



actually stabilizes the Kondo spin liquid phase, resulting 
in monotonic increase of both the spin and charge gaps 
as (X function of U c . 

The model we consider in this article is the fol- 
lowing one-dimensional KLM with Coulomb interaction 
(KLMC): 

H = -*E(4^+i CT + h.c.) + j5»f 

ier ifj, 



where c- 



{ci a ) is the creation (annihilation) operator 



of a conduction 
(1/2) 4^ 



electron at the ith site, and at 



' aa' 



x,y, z) are the spin density operators of the conduction 
electrons. In the second term represents a localized 
/-spin operator with S = 1/2. Concerning the kinetic 
energy of the conduction electrons we consider only the 
hopping processes between the nearest-neighbor sites, t. 
This model was considered by Yanagisawa and Harigaya 
in connection with the ferromagnetic ground state in the 
strong-coupling limit of the one-dimensional KLM away 
from half-filling B . Here we are interested in the half fill- 
ing case where the total number of conduction electrons is 
equal to the number of lattice site: N = ^2 ia . Ci a Ci a = L. 
It is needless to say that this Hamiltonian is reduced to 
the Hubbard model in the limit of vanishing exchange 
interaction, J — > 0, and to the usual KLM for U c = 0. 

For the case of U c = extensive studies on the spin and 
charge gaps have already been done. First Tsunetsugu 
et al. showed numerically that the ground state is singlet 
and the spin gap always exists for any finite exchange 
by usingexact diagonalization combined with finite-size 
scaling ||. They have concluded that the J dependence 
of the spin gap has a similar form as the Kondo temper- 
ature of the single impurity Kondo model but the expo- 
nent is larger than that of single impurity case, the lattice 
enhancement effect. Subsequently, semiclassical analysis 
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based on the nonlinear a model obtained by mapping 
from the one-dimensional KLM was carried out by Tsve- 
lik ||. He showed that the spin gap exists for any J in 
consistent with the above results with some logarithmic 
correction in the exponent. For the charge gap, on the 
other hand, a J-linear behavior is suggested by Nishino 
and Ueda. By using exact diagonalization of the one- 
dimensional periodic Anderson model up to eight sites, 
they observed that the charge gap is proportional to 1/U 
in the region of strong Coulomb interaction which means 
J-linear dependence in the corresponding weak-coupling 
region of the KLM (J = 8V 2 /U) §. It follows that the 
ratio between the charge and spin gaps A c /A s diverges in 
the limit of vanishing exchange constant. Yu and White 
applied the DMRG method to the KLM @. Their results 
for larger clusters have confirmed the diverging behavior 
of A c /A,5 in the region of weak exchange coupling. For 
an asymmetric case, however it is also pointed out the 
ratio remains finite Q . 

In spite of these progresses, the results in the region of 
weak exchange coupling are still not accurate enough to 
determine precisely the functional form of the gaps due to 
the smallness of clusters used in these numerical studies. 
For example, the exponent of the spin gap has not yet 
been determined with sufficient accuracy and the linear 
J dependence of the charge gap remains to be a plausible 
scenario. In order to determine functional forms of these 
gaps it is essential to treat large clusters systematically 
so that we can extrapolate the gap energies to the bulk 
limit. For this purpose the DMRG method developed by 
White J^] is efficient because the number of the states 
used to construct the wave function does not increase 
with increasing the system size. Naturally the trunca- 
tion of the states introduces numerical errors, but the 
errors may be estimated by the eigenvalues of the den- 
sity matrix which are truncated off. Thus it is possible 
to increase the system size within a given accuracy. 

In the previous DMRG calculation by Yu and White 
|^| , the number of states kept for each block is 180 and 
the number of the maximum system size is 24. However 
the lattice sizes are too small to fix the J dependence of 
the spin and charge gaps. For this purpose we use the 
finite system algorithm of the DMRG method with open 
boundary conditions keeping up to 300 states for each 
block and increase the system size up to 80. 

In contrast to the infinite system algorithm, the finite 
system algorithm gives more accurate results, but the 
extrapolation to the infinite system is necessary. The 
extrapolation to the infinite may be done in the following 
way. Since the lowest excited state generally corresponds 
to the bottom of an excitation spectrum which can be 
expanded in terms of k 2 , we expect finite-size scaling is 
of the form of L~ 2 for large systems. It means the gap 
energy of the large systems behaves as 

A(L) = A(oo)+ (3L- 2 + 0(L- 4 ). (2) 

We determine the gap energy of the bulk system by using 
this scaling when data for available system sizes already 



follows the scaling form. When a gap is tiny, for example, 
the spin gap of J = 0.6t and U c = 0, the scaling form is 
not yet clear for the system size available. In such a case 
the upper and lower bounds are estimated by L~ 2 and 
L _1 scaling, respectively, because the size dependence of 
the gap for smaller systems is close to L _1 rather than 
L- 2 . 

First we discuss the spin gap. It is obtained from the 
difference of the ground-state energies in the subspace of 
total S z being zero and one with the same total electron 
number L, E g (S z = 1,N = L)-E g {S z =0,N = L). The 
SU(2) symmetry in the spin space guarantees the energy 
difference is the same as the spin gap in the subspace of 
total S z being zero. 

Before proceeding to the analysis of the results it is 
worth reminding the Kondo impurity model. In the 
Kondo impurity model all physical quantities are scaled 
by single energy scale Tk ~ D exp which is known 

as the Kondo temperature. Here p — is the density 
of states of the conduction band at the Fermi level. The 
simplest extension to the lattice problem is to include an 
enhancement of the exponent owing to intersite correla- 
tions. Then the spin gap is expected to behave as 

A s ocexp( — ), (3) 

apJ 

where the a is the enhancement factor. The Gutzwillcr 
approximation predicts an enhancement factor of 2 ||. 
As a matter of fact, Tsunetsugu et al. have estimated 
that the enhancement factor a is in the range of 1 < 
a < 5/4 by using a finite-size scaling based on exact 
diagonalizations for the systems up to L = 10. 

The first task of the present study is to determine this 
enhancement factor more precisely. For this purpose we 
plot logarithm of the spin gaps as a function of 1/ J. Fig- 
ure [l] shows the results extrapolated to the bulk limit 
using data of L = 6, 8, 12, 18, 24, 40. From this figure we 
obtain the exponent a = 1-4(1) for the case of U c = 0. 
There are some ambiguities in the extrapolation to the 
bulk limit for a tiny gap. However within the present 
accuracy we do not observe any logarithmic correction 
to the exponent which is predicted by the semiclassical 
approach. 

Now we proceed to the effect of the Coulomb interac- 
tion between conduction electrons. As is well known, one 
of the most important effect of the Coulomb interaction 
for low-energy properties is the mass enhancement, which 
may be represented by an effective density of states at the 
Fermi energy, p* . Thus in the weak-coupling region it is 
natural to extend the form of Eq. (^) to finite U c by al- 
lowing U c dependence of the exponent, a(U c ). In Fig. |l| 
it is seen that the numerical data are nicely fitted by this 
form. The exponent a is determined for various U c . As 
shown in Fig. |^ the exponent a increases with increasing 
U c and the asymptotic behavior is linear in U c . 

In the limit of strong Coulomb interaction (U c /t — > oo) 
the KLMC is mapped to the Heiscnbcrg chain coupled 
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with localized / spins. Since the effective coupling of the 
Heisenberg chain is given by J c g = At 2 /U c , we can com- 
pare the present spin gap with that of the spin system. 
From Fig. |^ we find the asymptotic behavior of the expo- 
nent is a = 0.78U c /t + 0.7. It means that the spin gap of 
the above spin system behaves as A s ~ exp (— 2 J c g/ J). 
In order to check this form we have analyzed the numer- 
ical data for the spin system obtained by Igarashi et al. 
p0{ and found good agreement. Thus we conclude that 
the spin gap of the KLMC always vanishes exponentially 
in the limit of weak exchange coupling for any Coulomb 
interaction, U c , and the exponent a increases monotoni- 
cally with increasing U c . 

In the strong-coupling region of J, we may use the per- 
turbation theory with respect to tj J also in the presence 
of finite Coulomb interaction U c . After straightforward 
calculation we get the following result for the spin gap: 



A, = J 



At 2 



2t 2 



±J + U C p + u c 



(4) 



The solid curve in Fig. [I] represents this result for U c — 
lOi which shows that our numerical data for the spin 
gap is well reproduced by the perturbation result down 
to J ~ t. From Eq. (Q) one can see that the deriva- 
tive of the spin gap with respect to U c is always positive, 
A s > 0, for any positive U c and J. From these obser- 
vations in both the strong- and weak-coupling regions we 
can conclude that the spin gap increases with increasing 
Coulomb interaction for any antiferromagnetic exchange 
J. 

The charge gap is obtained by E g (S z = 0, N = L + 
2) - E g (S z = 0,N = L). Owing to the hidden SU(2) 
symmetry in the charge space, the energy difference is 
the same as the charge excitation gap in the subspace of 
total electron number fixed to L 0. 

Before discussing the charge gap we first notice the 
relation between the charge gap A c and the quasipar- 
ticle gap A qp which is defined by E g (S z = ±1/2, N — 
L ± 1) - E g (S z = 0,N = L). In the strong-coupling 
limit, J /t — > oo, it is evident that the charge gap is twice 
the quasiparticle gap because the energy required to cre- 
ate the lowest charge excited state is the same as the 
energy cost to add two additional electrons owing to the 
SU(2) symmetry in the charge space. In the second order 
perturbation in t/J, one can show that the interaction 
between the two additional electrons is repulsive, leading 
to only a phase shift. Therefore, the charge gap is given 
by the sum of two quasiparticle gap A qp in the bulk limit: 



A r = 2A 



qii- 



(5) 



A similar argument is also valid for the periodic Ander- 
son model H . Validity of this relation is checked by the 
present DMRG calculation and we have confirmed this 
relation in the entire region of the exchange constant J. 
On the other hand the spin gap is determined by the 
lowest bound state of a quasielectron and a quasiholc. 



Let us start from the case of J = where exact results 
are known. In this limit the KLMC is reduced to the 
Hubbard model which is exactly solved by Lieb and Wu 
for the one-dimensional cas e The asymptotic form 
of the charge gap is A c cx y/U c t exp ( — jj-) for small U c , 
and A c cx U c — At for large U c . Figure [jjshows the charge 
gap extrapolated to the infinite system from the data 
for L = 6,8, 12, 18, 24, 40. These results for finite J are 
consistent with the exact results which are denoted by 
crosses on the vertical axis. 

For U c — 0, the charge gap of the Hubbard model 
vanishes. In this case an important question is how the 
charge gap opens for finite J. As is shown in the previous 
work by Nishino et al. the charge gap is much bigger 
than the spin gap in the weak-coupling regime j5|. It 
implies that the correlation length for the spin degrees 
of freedom is much longer than the charge correlation 
length. Therefore, for the discussion of the charge gap 
it is justified to assume that the spin-spin correlation 
length is infinitely long. Under the assumption of the 
infinite spin correlation length, one can get the charge 
gap which is linear in J with its coefficient 1/2 



J 
2' 



(6) 



This linear dependence in J is actually seen in the present 
DMRG calculations shown in Fig. For the case of 
J = 0.2t and U c = we have additionally calculated a 
larger system of L = 80 and we get the charge gap of 
O.lt from the extrapolation. We may take this fact as a 
numerical confirmation of the coefficient of 1/2. It should 
be stressed again that this type of mean-field theory is 
not justified for the discussion of the spin gap. 

From Fig. || we also find that the charge gap increases 
with increasing Coulomb interaction U c . In the strong- 
coupling regime, we get the following result for the charge 
excitation gap within the second order perturbation: 



-J + U c -2t 



-A 2 



IJ + U C 



3*2 
J ' 



(7) 



From the result it is clear that the charge gap is a in- 
creasing function of U c for large J. Thus it is concluded 
that the charge gap increases with increasing Coulomb 
interaction U c for any exchange constant J similarly as 
the spin gap. 
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FIG. 1. Spin gap of the one-dimensional Kondo lattice 
model with Coulomb interaction. The thick curve repre- 
sents the result of the perturbation theory in terms of t/J 
for I7 C = Wt. Typical truncation error in the DMRG calcu- 
lation is 10 -6 for J = 1. Errorbars are estimated from L^ 1 
and L~ 2 scaling. Gap energies, exchange constant J, and 
Coulomb interaction U c are in units of t. 

FIG. 2. (7 c -dependence of the exponent of the spin gap. 
Coulomb interaction U c is in units of t. 

FIG. 3. Charge gap of the one-dimensional Kondo lattice 
model with Coulomb interaction. Results on the vertical axis 
are obtained from the exact solution of Lieb-Wu. Typical 
truncation error in the DMRG calculation is 10 -6 for J — 1 
and 10 - for J = 0.2, which is a dominant source of numerical 
errors since the finite-size scaling, Eq. (2), is well obeyed. Gap 
energies, exchange constant J, and Coulomb interaction U c 
are in units of t. 
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